en 
(N 



THE FRACTIONAL POISSON PROCESS AND THE INVERSE 

STABLE SUBORDINATOR 

^H ! MARK M. MEERSCHAERT, ERKAN NANE, AND P. VELLAISAMY 

O; 

Abstract. The fractional Poisson process is a renewal process with Mittag-Lefher 

""j" , waiting times. Its distributions solve a time-fractional analogue of the Kolmogorov 

fT ^ ' forward equation for a Poisson process. This paper shows that a traditional Poisson 

process, with the time variable replaced by an independent inverse stable subordina- 

tor, is also a fractional Poisson process. This result unifies the two main approaches 

in the stochastic theory of time-fractional diffusion equations. The equivalence 

I— I' extends to a broad class of renewal processes that include models for tempered 

Ph , fractional diffusion, and distributed-order (e.g., ultraslow) fractional diffusion. The 

Mh ' paper also establishes an interesting connection between the fractional Poisson pro- 

(~| , cess and Brownian time. 

-)— > . 
C^ . 

B: 

' ^; 1. Introduction 

^ ■ The fractional Poisson process (FPP) was introduced and studied by Repin and 

^~~^ '. Saichev |12], Jumarie [22], Laskin [27], Mainardi et al. [301 [32], Uchaikin et al. [17] 

Q I and Beghin and Orsingher [5l E]. The FPP is a natural generalization of the usual 

lO ' Poisson process, with an interesting connection to fractional calculus. This renewal 

t~>^ ■ process has IID waiting times J„ that satisfy 

O: (1-1) F{Jn>t)=E^{-Xt^) 

. . ! for < B < 1, where 

> . 

• w^ °° k 



fc=0 

denotes the Mittag-Leffler function. When (3 = 1, the waiting times are exponential 
with rate A, since e^ = Ei{z). Let T„ = Ji -|- ■ ■ ■ -|- J„ be the time of the nth jump. 
Then the FPP 

(1.3) Np{t) = max{n > : T„ < t} 

is a renewal process with Mittag-Leffler waiting times. 

A compound FPP is obtained by subordinating a random walk to the FPP. The 
resulting process is non-Markovian (unless /3 = 1) and the distribution of that process 
solves a "master equation" analogous to the Kolmogorov equation for Markov pro- 
cesses, with the usual integer order time derivative replaced by a fractional derivative. 
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The continuous time random walk (CTRW) is another useful model in fractional 
calculus. Consider a CTRW whose IID particle jumps Yn have PDF w{x), and whose 
IID waiting times (J„) are Mittag-Leffler variables independent of (1^„). The particle 
location after n jumps is S{n) = Yi + ■ ■ ■ + Yn, and the CTRW S{Ni3(t)) gives the 
particle location at time t > 0. Hilfer and Anton [19] show that the PDF p{x, t) of 
the CTRW S{Ni3{t)) solves the fractional master equation 

/oo 
p{x-y,t)w{y)dy 
-oo 

where 9f denotes the Caputo fractional derivative. The Caputo fractional derivative, 
defined for < n — 1 < /3 < n by 

(1.5) d^g{t) = ^^J_^^ I* (t - rr-'-^9^''\r)dr, 

where g^''^ denotes the k-th derivative of g, was invented to properly handle initial 
values [TT] . 

If /3 = 1, then (9f is the usual first derivative. The corresponding CTRW S{Ni{t)) 
is a compound Poisson process, and ( 11. 4p reduces to 

/oo 
p{x-y,t)w{y)dy, 
-oo 

the Cauchy problem associated with this infinitely divisible Levy process. Then a 
general result on Cauchy problems [2], Theorem 3.1] implies that the PDF of the 
time-changed process S{Ni{E(t))) solves the fractional Cauchy problem (11.41) . where 

(1.7) E{t) = inf{r > : D{r) > t} 

is the right-continuous inverse (hitting time, first passage time) of D{t), a standard 
/3-stable subordinator with E[e~'^^'^*''] = e"*" for some < /3 < 1. 

Since the PDF of both 5'(iVg(t)) and S{Ni{E{t))) solve the same governing equa- 
tion (11.41) . with the same point-source initial condition (i.e., both processes start at 
the origin), these two processes have the same one dimensional distributions. Heuris- 
tically, the degenerate case Yn = I gives S{n) = n, which strongly suggest that the 
FPP iV^(t) and the process Ni{E{t)) have the same one dimensional distributions. 
We will call Ni{E{t)) the fractal time Poisson process (FTPP), since it comes from a 
self-similar time change (see, e.g., [SU Proposition 3.1]). In this paper, we will prove 
that the FPP and the FTPP are in fact the same process, by showing that the wait- 
ing times between jumps in the FTPP are IID Mittag-LefHer. This strong connection 
between the FPP and the FTPP unifies the two main approaches in the stochastic 
theory of fractional diffusion. For example, the FPP approach was used recently in 
the work of Behgin and Orsingher [5] , while the inverse stable subordinator is a key 
ingredient in |38j . 
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2. Two EQUIVALENT FORMULATIONS 

Recall that the fractional Poisson process (FPP) Nf^{t) is a renewal process with 
Mittag-Leffler waiting times fll.ip . and the fractal time Poisson process (FTPP) 
Ni{E{t)) is Poisson process, with rate A > 0, time-changed via the inverse stable 
subordinator (11. 7p . The proof that the FPP and the FTPP are the same process 
requires the following simple lemma. 

Lemma 2.1. Let D{t) he a strictly increasing right- continuous process with left-hand 
limits, and let E{t) be its right- continuous inverse defined by fll.7p . Then 

(2.1) D{r-) = sup{t > : E{t) < r} 
for any r > 0. 

Proof. Let to = sup{t > : E(t) < r}. Then there exists a sequence of points t„ f ^o 
such that E{tn) < r for all n. Let e„ = r — E{t„,) > 0. If r > E{t) then, since 
D{t) is strictly increasing, D{r) > t. Since D{r) is right-continuous, it follows that 
L'(E(t)) > t for all t > 0. Then we have t„ < D{E{tn)) = D{r - En) < D{r-). 
Letting n — )■ oo shows that D{r—) > to- 

Since D has left-hand limits, for any r„ f r we have D{rn) — )■ D{r—) as n — )■ oo. If 
D{r—) > to, then for some rn < r we have D{rn) > to- Since E{t) is nondecreasing 
and continuous, this implies that E{D{rn)) > r, by definition of to. But, E{D{r)) = r 
for all r > implying that r^i > r, which is a contradiction. Thus, (12. ip follows. D 

Theorem 2.2. For any < (3 < 1, the FTPP Ni{E{t)) is also a FPP. That is, the 
waiting times between jumps of the FTPP are IID Mittag-Leffler. 

Proof. Let W„ be an IID sequence with P(W„ > t) = e"-^* and K = W^i H VWn 

so that the Poisson process Ni{t) = maxjn > : \4 < t}. Let 

(2.2) Tn = sup{t > : Ni{E{t)) < n} 

denote the jump times of the FTPP. This definition of the jump times takes into 
account the fact that E(t) has constant intervals corresponding to the jumps of the 
process D{t). Using the fact that {Ni{t) < n} = {Vn > t} for the Poisson process, 
along with (12. 2p . we have 

r„ = sup{t > : E{t) < K}- 

Then Lemma l2.ll implies that Tn = DiVn—). Define Xi = ti and Xn = Tn — Tn-i 
for n > 2, the waiting times between jumps of the FTPP. In order to show that the 
FTPP is an FPP, it suffices to show that Xn are IID Mittag-Leffler, i.e., they are IID 
with Jn- 

Recall that the Laplace transform of the exponential distribution E(e~*^") = 
A/(A -I- s). Also recall that E(e~'*'^*^*)) = e~** . Since D{t) is a Levy process, it 
has no fixed points of discontinuity and hence D{t—), D{t) are identically distributed 
for all t > 0. (Indeed, D{t) = D{t~) a.s. [H Lemma 2.3.2]). 
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(2.3) 



Then a conditioning argument yields 



-Wisf 



X 



A + s/3 



Let ff}{x) = d^[l - Ep{-Xx^)] be the Mittag-Leffler PDF of J„. It is well known that 



'Esi-Xx'^) dx 



-3-1 



'P 



X + s^' 
see for example [38| Eq. (3.4)]. Now integrate by parts to see that 

fOO 



(2.4) 



E(e-'^^) = / e-'''fp{x)dx 
Jo 

= / se""^ (1 - Ei3{-Xx^)) dx 
Jo 



,/3-l 



A + s/^ 



A 



X + s^ 



E(e" 



and then the uniqueness theorem for LT implies that Ti, ri are identically distributed. 
In particular, Xi has the same Mittag-Leffler distribution as Ji. 

A straightforward extension of this argument shows that (Ti, . . . , T„) is identically 
distributed with (ri,...,r„) for any positive integer n. To ease notation, we only 
write the case n = 2. First observe that 

E(e--i^ie-^^^^) = E(e-^^'^^e-*^('^^+-^^)) = -^ -^ ■ ^--,, 

A+(Si + S2)^ X + s'2 

using the independence of Ji and J2. Next write 

= E(e"(''i+'*2)^(*i)e"*2[^(*^+*2)"^(*^)^) 

using the fact that D{t) has independent increments. Then 

]g('g-SlTi-S2T2\ _ ]g/'g-SlZ)(Vl/i-)-S2£'{[Wi+Ty2]-)\ 



E 



-H/i(si+52)'5^-W/2sf 



A 



A 



E(e-''^'e-'^^'] 



X + {si + S2y A + 4 

Now an application of the continuous mapping theorem shows that (Ji, . . . , J„) is 
identically distributed with {Xi, . . . ,Xn) for any positive integer n. Then (X„) is an 
IID sequence, so Ni{E{t)) is a renewal process. D 
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Remark 2.3. Theorem 12.21 extends a result in Behgin and Orsingher [S]. They define 
(in our notation) a random variable E{t) and show that the two random variables 
Ni^it) and Ni{E{t)) have the same density function, by comparing their Laplace 
transforms. They identify E{t) only though its density function, which they express 
in terms of an integral involving the density of D{t), see Remark 13.21 for more detail. 
Cahoy, Uchaikin, Woyczynski [10] also connect the Mittag-LefHer distribution with a 
stable law. They note that (in our notation) P{Jn > t) = E[exp{—Xt^/D{l)'^)], which 
is useful in simulations. To connect this with our work, note that E{t) = (t/D{l))^ 
in distribution (see Corollary 3.1 in [M]), so that P(J„ > ^) = E[exp{—\E{t))]. A 
result of Bingham [S] shows that the the Laplace transform of the stable hitting time 
E{t) is Mittag-Leffler, so that (dH]) holds. 

Remark 2.4. The proof of Theorem 12.21 uses the fact that, if -D(t) is a /3-stable subordi- 
nator and Wi is exponential, then D{Wi) has a Mittag-Leffler distribution. This fact 
was first noticed by Pillai [ID], who showed that W^ D(l) is Mittag-Leffler. These 
are equivalent because D{t) is identically distributed with t^^^D{l). This Mittag- 
Leffler distribution is also known as the positive Linnik law, e.g., see Huillet [21]. It 
has the property of geometric stability: A geometric random sum of Mittag-Leffler 
variables is again Mittag-Leffler, e.g., see Kozubowski 



Next we want to show that the FTPP Ni{E{t)), and hence also the FPP iV/3(t), 
occurs naturally as a CTRW scaling limit. This provides a further justification for 
the FPP as a robust physical model, see for example Laskin [27]. Suppose now that 
P(Jn > t) = t~^L(t), where < /3 < 1 and L is slowly varying. For example, this 
is true of the Mittag-Leffler waiting times. Then Ji belongs to the strict domain of 
attraction of some stable law D with index < /3 < 1, i.e., there exist &n > such 
that 

(2.5) 6„(Ji + --- + J„,)^D, 

where D{1) = D > {] almost surely, and ^ denotes convergence in distribution. Let 
h{t) = b[t]. Then b{t) = t~^/^LQ{t) for some slowly varying function Lo(t) (e.g., see 
[T71 XVII. 5]). Since h varies regularly with index — 1//3, h~^ is regularly varying with 
index 1//3 > and so by [^ Property 1.5.5] there exists a regularly varying function 
h with index /3 such that l/6(6(c)) ~ c, as c — )■ oo. Here we use the notation f ^ g for 
positive functions /, g if and only if f{c)/g{c) — )■ 1 as c ^ oo. Let T„ = Ji -(- • — \- Jn 
and define a renewal process 

(2.6) R{t) = max{n > : T„ < t} 

with these waiting times. Next, construct a CTRW with iid Bernoulli jumps Yn 
with P(ri^^ = 1) =P and P(yi^^ = 0) = 1 - p, independent of (J„). Let S^P^n) = 
y/^ + ■ ■ ■ + Yn , a binomial random variable. Then S^P\R{t)) is a CTRW with heavy 
tailed waiting times and Bernoulli jumps. 
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Theorem 2.5. The FTPP is the process limit of a CTRW sequence: 
(2.7) {S^"'^''^\[\R{ct)])]^^^ => {N,{E{t))}^^^ 

as c —)■ oo in the Mi topology on D{[0, oo),M). 

Proof. Since the sequence (J„) is in the strict domain of attraction of a /3-stable 
random variable D, [33] CoroUary 3.4] shows that 

{b{c)-'R{ct)}^^^ => {E{t)}^^^ as c ^ oo. 

in the Skorokhod Ji topology, where D{t) is the stable subordinator with -D(l) = D, 
and E{t) is given by (11. 7p . 

Since the binomial random variable S^^\n) has LT E(e~'^'^ '' ^^^) = (1 + (e^'^ — l)p)"' 
for any n > 0, it follows that 

^^g-.5W([At/p])) ^ (1 ^ (g~. _ I)p)[^*/P] ^ exp(-At(l - e-^)), 

as p — )■ 0, using the fact that (1 + apY^^ — )■ e*^ as p — )■ 0. It follows by the continuity 
theorem for LT that S^^\[\t/p]) =^ Ni{t) for any t > 0, since exp(— At(l — e~*)) is 
the LT of the Poisson random variable Ni{t). Then a standard argument (e.g., see 
[33| Example 11.2.18] shows that we also get 

{S''\[^t/p])h>o^{Ni{t)},>o, 

as p -T- 0, where ==^ denotes convergence of all finite dimensional distributions. Since 
the sample paths of S^'P\[\t/p]) are increasing and Ni{t) is continuous in probability, 
being a Levy process, Ji convergence follows using [SI Theorem 3]. 

Since the CTRW waiting times ( J„) are independent of the jumps {Yn ), and since 
l/6(c) — ;• as c — !■ oo, it follows that 

{S^'''Kc))^^Xth{c)]):h{c)-'R{ct)) => (iVi(t),i?(t)) 

in the Ji topology of the product space -D([0, oo), M x M), by [3 Theorem 3.2]. Since 
the process E{t) is nondecreasing and continuous, [15| Theorem 13.2.4] along with 
the continuous mapping theorem yields 

SW''('^^\[XR{ct)]) = S(^/''('^))([A ■ b{c)-^R{ct) ■ b{c)]) => Ni{E{t)) 

in the Mi topology on D([0, oo), M). D 

Remark 2.6. For the specific case of Mittag-Leffler waiting times, where P( Jn > ^) = 
E/si—t^), we can take 6„ = n~^^^ in (12. 5p . To check this, note that 
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as n — 7- oo. Then 6(c) = c^ and the CTRW convergence (12.71) reduces to S^'^ \[\R{ct)\) 
Ni[E{t)) as c — )■ oo. Substitute p = c~^ to get 

(2.8) S^P^lXRip-^^H)]) ^ Ni{E{t)), asp^O. 

3. Fractional calculus 

This section develops some interesting connections between the fractional Poisson 
process and fractional calculus. In the process, some apparent inconsistencies in the 
existing literature will be explained. Behgin and Orsingher Eq. (2.17)] show that 
the FPP of order < /3 < 1 has distribution 



(Ax)'= 



(3.1) W>(NB(t) = k)= e-^''^—^V(x,t)dx, 

Jo k\ 

where V{x, t) is a "folded PDF" defined on a; > 0, for each t > 0, by V{x, t) = 2v{x, t), 
and v{x, t) is another PDF with x G M for each t > that solves 

dfv{x,t) = dlv{x,t); 

(3.2) v{x,0) = 6{xy, 

dtv{x,0) = 0, a 1/2 < (3 < I. 

It is also stated in [51 Eq. (1.9)] that the FPP iV/3(t) = Ni{Tt), where Tt is a random 
process with PDF V{x,t) for t > 0. However, that process is identified only in terms 
of its one dimensional distributions (PDF). Theorem 12.21 shows that the inverse stable 
subordinator E{t) is one such process. 

On the other hand, a simple conditioning argument shows that the equivalent 
FTPP process has distribution 

/•oo POO (Xr)'' 

(3.3) F(Ni(E(t)) = k) = P(Ni(x) = k)h(x,t)dx= e-^''^-^h(x,t) dx 

Jo Jo k\ 

where h{x,t) is the density of E{t), a PDF on x > for each t > 0. It follows from 
[361 Theorem 4.1] that this PDF solves 



(3.4) d^h{x, t) = -d^h{x, t); h{x, 0) = 6{x). 

In view of Theorem 12.21 the two distributions (13. ip and (13. 3p must be equal. Thus, 
the main purpose of this section is to reconcile the two fractional differential equations 
(13:21) and 1^. 

Theorem 3.1. Let Nis{t) be a fractional Poisson process (11.31) with < /3 < 1, so 
that (13. ip holds. Let Ni{E{t)) be the equivalent fractal time Poisson process, where 
E{t) is the standard inverse (3-stable subordinator with PDF h{x,t), so that (13.31) 
holds. Then 

(3.5) h{x,t) = 2v{x,t) for all X > andt> 0. 
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In particular, the two fractional partial differential equations fl3.2p and fl3.4|] are con- 
sistent, in the sense that the folded solution V{x,t) = 2v{x,t) to (I3.2p coincides with 
the solution h{x,t) to (13. 4p . 

Proof. Mainardi [29l Eq. (3.2)] shows that the solution to the fractional diffusion-wave 
equation (13. 2p has LT 

(3.6) v{x,s)= e~''v{x,t)dt = -s^^^e~\^\'^ 

Jo 2 

while [Ml Eq- (3.13)] shows that 

(3.7) h{x,s) = s^-^e-'^''. 

Since both are different iable in t, they are also continuous, so LT uniqueness for 
continuous functions implies (13. 5p . 

Take Fourier transforms in (13. 6p to see that the solution to (13. 2p has Fourier-Laplace 
transform (FLT) 

roc roc £;2/3— 1 

(3.8) v{k,s)= e-'* e-'^''v{x,t)dxdt 



s'^ + k'' 
where we have used the fact that e~"'^' has FT 2a/(a^ + /c^). Rearrange to get 

s^^v{k,s)-s^^-' = -ev{k,s) 
and invert the FT to get 



s'^^vix, s) — s 



2^-^^;(x,0) =dlv{x,s), 



using the fact that dxf{x) has FT {ik)f{k) and f(x,0) = 5{x) has FT v{k,Q) = 1. 
To invert the LT, note that the Caputo fractional derivative d^ f{t) has LT s^f(s) — 
s^-7(0) if < /3 < 1, and LT s^f{s) - s^-^f{0) - s^-^f'{0) if 1 < /3 < 2. This 
is easy to verify from the definition (II. 5p . using the corresponding formula for the 
integer derivative, along with the fact that s^~^ is the LT of t~''/r(l — /3). Now use 
the remaining initial condition dtv{x,0) = for 1/2 < /3 < 1 to invert the LT, and 
arrive at (13. 2p . 

Likewise, the solution to (13. 4p has FLT 

(3.9) h{k,s)= e-'^"/-^e-"^^dx 



JO 31^ + ik' 

using the fact that e°'^I{x > 0) has FT l/(a + ik). To see that these are consistent, 
compute the FLT of h{\x\,t): 

/OO /"OO / /"OO /»oo \ 

e-'^''h{\x\,t)dxdt= e"'M/ e-'''^h{x,t)dx+ e'^''h{x,t)dx] dt 

-OO Jo \Jo Jo J 

«/3-l _g/3-l / 52/3-1 \ 



sP + ik sP-ik \s'^f^ + k^ 



Invert the FLT to see that h{\x\,t) = 2v{x,t) for all a; G M and t > 0. To verify the 
LFT solution, take FT in (13. 4p to get 

dfh{k,t) = -ikh{x,t) 

and apply the LT to get s^h{k^ s) — s^~^ = —ik h{k, s), using the point source initial 
condition h{k,0) = 1. D 

Remark 3.2. Behgin and Orsingher [5l Eq. (2.21)] show that 

v{x,t) = 2Y(l-R) J ^^"^) ^p{\xlt)dx 

where p{x,t) is the density of the stable subordinator D(t), while [361 Theorem 3.1] 
implies that 

1 /"* 

h{x,t) = Y(i_ R) / *^^^^) ^p{x,t)dx. 

This gives another proof that h{\x\,t) = 2v{x,t). 

Remark 3.3. A closed form expression for the probability mass function p{n,t) = 
F{Ni{E{t)) = n) = P{Np{t) = n) follows easily from (Q. Use |36l Eq. (3.13)] to 
write 

"' dx 





p{n, s)= e '^p{n, 


t)dt= 1 


1^- 


and us€ 


; the formula for the gamma 


density 


to compute 


(3.10) 


p{n,s) 


5/3-] 


A" 


A + .' 


i/3(A + s/5)"' 


Invert i 


ising the generahzed Mittag 


-Leffler function 




ElA^ 


00 

r=0 


(7)r^^ 

r\T{ar + 9) 



where (7)^ = 7(7 + 1) ■ ■ ■ (7 + r — 1) is the Pochammer Symbol. Formula (2.5) of 
gives 

/•CO a'^S-'J 

Substitute z/ = /3, (5 = n + l and 7 = /3n + 1 to get 

(3.11) p(n.O ^ (M^TErU^-Xt^) J^t^^-±^y 

which is the same form obtained by Jumarie [22], Laskin [27], Beghin and Orsingher 
[3 [6] and Cahoy [9] using different methods. 
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The equivalence in Theorem 13.11 results from folding the solution to the fractional 
diffusion- wave equation fl3.2p . Another fractional partial differential equation for the 
density h{x,t) of the standard inverse /3-stable subordinator E(t), which is closer to 
the form (13. 2p . can be obtained by arguments similar to those used in [3] to connect 
the inverse stable subordinator to iterated Brownian motion. In that theory, it is 
customary to avoid distributions by imposing a functional initial condition. 

Theorem 3.4. Let E{t) be the standard inverse P-stable subordinator with density 
h{x,t). Then for any f G L2(M) fl C^(M), the function 

POO 

(3.12) u{x,t)=E^[f{E{t))]= / fix + y)h{y,t)dy 

Jo 

solves the fractional differential equation 

(3.13) dfu{x,t) = -dj{x) ^^^^^^ +dlu{x,ty, u{x,0) = fix). 
In particular, when /3 = 1/2, (13.121) solves 

(3.14) dtu{x,t) = ~^''f}^> +dlu{x,ty, u{0,x) = f{x), 

VTTt 



and in this case we also haveu{x,t) = E,x[f{\B{t)\)], where B{t) is a Brownian motion 
with variance 2t. 

Proof. From (13.91) . we have 

- ^^"'/(^) _ g^~V>) S" -Zk _ g2/3-l _ ,fcg/3-l 
""^ '^^~ S^ + tk ~ sP+ik ' S^-lk~ s^f' + k^ ^^' 

SO that s^^u{k,s)-s'^f^-^f{k) = -ikf {k)s^~^ -k^u{k,s), which inverts to IKUh . It is 
well known that the Brownian motion first passage time D{y) = inf {t > : B{t) > y} 
is a stable subordinator with index /3 = 1/2 [H Example 1.3.19]. Then it is easy to 
see that 

E{t) = ini{y > : D{y) > t} = sup{5(r) : < r < t} 

and this recovers the fact, typically proven using the reflection principle, that 

P{E{t)>y) = 2PiB{t)>y). 

Then E(t) and |-B(t)| have the same one dimensional distributions, so we also have 
u{x,t) = Kx[f{\B(t)\)]. Note that X{t) = X{0) —t is a fortiori a continuous Markov 
process associated with the shift semigroup T(t)f{x) = Ej.[/(X(t))] = /(x — t) with 
generator 

i->-0+ t 
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Then [51 Corollary 3.4] implies that u{x,t) solves the equation 

dtu{x, t) = "^ + LJu{x, t); m(0, x) = f{x). 

Vvrt 



When L^ = —d^, this reduces to fl3.14p . a special case of flS.lSp with r(l/2) = ^/n. D 

Remark 3.5. In the case f{x) = S{x), Theorem 13.41 gives an alternative governing 
equation for h{x,t). Note that fl3.13p is very similar to the governing equation (13. 2p 
for the unfolded PDF. 

Remark 3.6. The process |-B(t)| in Theorem 13.41 is not the same process as the in- 
verse 1/2-stable subordinator E{t) in Theorem 12.21 although they have the same one 
dimensional distributions. Hence, the FTPP is not the same as the Brownian time 
Poisson process Ni{\B(t)\). However, we do have E{t) = sup{i?(r) : < r < t}, so 
that a Poisson process subordinated to the supremum of a Brownian motion is an 
FPP with (3 = 1/2. 

Remark 3.7. Let E{t) be the standard inverse stable subordinator of index /3 = 1/m 
for integer m > 1. Then [3l Remark 3.11], [391 Theorem 1.1] and Keyantuo and 
Lizama [23l Theorem 3.3] imply that u{x,t) = E,x[f{E{t))] solves 

'^-1 p/m.-l 

for t > and x G M, which is then equivalent to f l3.4p . The proof is similar to 
Theorem 13.41 For example, when /3 = 1/3 use 

""^^^ ^) - sl/3 + ik ' S^l^ - S^lHk + k^ ~ TTik^ ^^^' 

4. Renewal processes and inverse subordinators 

Theorem 12.21 shows that a Poisson process, time-changed by an inverse stable sub- 
ordinator, yields a renewal process with Mittag-Lefiler waiting times. This section 
extends that result to arbitrary subordinators that are strictly increasing. Let D{t) 
be a strictly increasing Levy process (subordinator) with E[e~^'^*-*^] = e"*'^^^*-', where 
the Laplace exponent 

Jo 
b > 0, and (po is the Levy measure of D. Then we must have either 

(4.2) 0d(O,oo) = oo, 

or 6 > 0, or both. Let E{t) be the inverse subordinator (11.70 . and recall that Ni{t) is 
a Poisson process with rate A. 
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Theorem 4.1. The time-changed Poisson process Ni{E{t)) is a renewal process 
whose IID waiting times (J„) satisfy 

(4.3) P(J„>t)=E[e-^^W]. 

Proof. The proof is similar to Theorem 12.21 Take Ni{t) = max{n > : \4 < t}, 
where Vn = Wi + ■ ■ ■ + Wn, with W„ HD as F{Wn > t) = e'^K Let 

r„ = sup{t > : A^i(^(t)) < n} = sup{t > : E{t) < K} 

and apply Lemma [2. II to get r„ = D{Vn—). Then, as in the proof of Theorem 12. 2 [ we 
have 



E(e~'^i) =E(e" 



^sD{Wi~)s 



^^•^) = E [E (e-^'^^^i^liyi)] = E [e-^i'^^W] ^ 



By [Ml Corollary 3.5], the IID random variables J„ in f l4.3p satisfy 

(4.5) l^-m, >t)M^ f .-'Ele-.'.].. ^ ^0^. 

Integrate by parts to get 

(4.6) ^ 

re-^'¥jM) = r se-^' [l-nJn >t)]dt = l- ^^I'l = ^ 

Jo Jo X + iPd{s) X + iIjd{s) 

which shows that Ti = Ji is identically distributed with ti. Extend this argument, as 
in the proof of Theorem 12. 2| to show that (Ti, . . . , T„) is identically distributed with 
(ti, . . . , Tn) for any positive integer n. For example, when n = 2, write 

and condition to get 

= E [E (e-^i^W)-^2^W+^2)|l^i,1^2)] 

A A 



E \q~^'>-'^o{si+S2) ^-W2i^D{s2)~\ 



A + ^d(Si + S2) A + V'd(s2) 

On the other hand, 

E(e-''^'e-''^') = E(e-''-^'e-''^-^'+-^'^) = -^ ■ ^— - 

X + iPd{si + S2) X + ijD{s2) 

using the fact that ( J„) are IID. To finish the proof, use continuous mapping to show 
that ( Ji, . . . , Jn) is identically distributed with {Xi, . . . , Xn), where Xn = Tn — t„_i 
are the waiting times between jumps for the process Ni{E(t)). D 

12 



Remark 4.2. Let A'^d(^) denote the renewal process from Theorem 14.11 so that 

(4.7) Noit) = max{n > : T„, < t}, 

where T„ = XlILi ^* ^^^ i-^n) are IID according to (14. 3p . Theorem 14.11 shows that 
Noit) = Ni{E{t)). This extends the relation A^^(t) = Ni{E{t)) from Theorem [Ml 
the special case of an inverse stable subordinator E(t) and Mittag-LefHer waiting 
times J„, to a general inverse subordinator. 

Remark 4.3. Let M{t) = E,{Nr){t)) denote the renewal function of the renewal process 
Nuit). Then using Lageras [261 Equation 4], it follows that the LT of M(t) is X/ipois). 

Remark 4.4. A simple conditioning argument shows that 

roo poo iXxY 

PD(n,t) = F(ND(t)=n)= PiNiix) = n)hix,t)dx = e'^^'^-^hix^t) dx 

Jo Jo nl 

where h{x,t) is the density of E{t). A straightforward extension of the argument in 
Remark 13.31 shows that 

which reduces to (I3.10p in the special case iPd{s) = s^ for a stable subordinator 
D{t). Use (USD and dHj) to see that the first factor in dHj) is the LT (t ^ s) of 
h{X, t) = F{Jn >t) = E(e-^^(*)), and the second is the LT of T„ = Ji + J2 + ••■ + Jn 
with Jrt IID as in (14. 3p . Denote the distribution of T„ by F^"''*\ the n-fold convolution 
of the distribution function F of Ji. Invert the LT to get 

(4.9) PDin,t)= f h{\,t-s)F^'''*\ds) 

Jo 

which extends (13. lip . 

5. CTRW SCALING LIMITS AND GOVERNING EQUATIONS 

In this section, we extend the fractional calculus results of Section [3] to the inverse 
subordinators of Section HI A general theory of CTRW scaling limits and governing 
equations is developed in [36]. Consider a sequence of CTRW indexed by a scale 
parameter c > 0. Take J^ nonnegative IID random variables representing the waiting 
times between particle jumps and T'^{n) = Yl^=i ^ii ^^^ time of the nth jump. Let 
Y^ be IID random vectors on W^ representing the particle jumps, independent of the 
waiting times, and set S^{n) = J2^=i ^i^ ^^^ location of the particle after n jumps. 
Define N^ = max{ra > : T'^{n) < t}, the number of jumps by time t > and 

(5.1) x\t) = s\N^) = y^y: 

i=l 
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the position of the particle at time t > and scale c > 0. Assume a triangular array 
limit 

(5.2) {{S'ict), r^(ct))},>o => {{A{t), D{t)}t>o, as c ^ oo, 

in the Ji topology on i5([0, oo),M'^ x IR+), so that A{t) and D(t) are independent 
Levy processes on W^ and M, respectively. Since the waiting times are nonnegative, 
D(t) is a subordinator. In this section, we assume the drift 6 = in fl4.ip . as well as 
condition f l4.2p and 

(5.3) / y\ In y\(f)D{dy) < oo. 

Assumption (14. 2 p implies that the process {D{t)} is strictly increasing, i.e., D{t) is 
not compound Poisson. Then [361 Theorem 3.1] shows that the inverse subordinator 
E(t) in (ll.7p has a Lebesgue density 

(5.4) h{x, t)= I 4>D{t - y, oo) FDi,){dy). 

Jo 

Write E[e-^-^W] = e-*'^^^'), as before. Let P{x,t) = F{A{t) < x) be the distribution 
function of A(t), and write 



p(k,t) = f e-'''''P{dx,t) =e" 



where ipA{k) is the Fourier symbol of A. The symbols define pseudo-differential 
operators: ipD{dt)f(t) has LT ipD{s)f{s), and ipA{—iDx)f{x) has FT ipA{k)f{k), for 
suitable functions /. Then |36l Theorem 2.1] establishes the CTRW scaling limit 

(5.5) {X'{t)}t>o => {A{E{t))}t>o, as c ^ oo, 

in the Mi-topology on D{[0, oo),M.'^). Recall that a function Q is a mild solution to 
a space-time pseudo-differential equation if its (Fourier-Laplace or Laplace-Laplace) 
transform solves the equivalent algebraic equation in transform space. The next result 
is a small extension of [361 Theorem 4.1]. 

Theorem 5.1. Assume (15. 2 p holds, where D{t) is a subordinator without drift such 
that conditions (14. 2 p and (15. 3p hold. The distribution function of the CTRW limit 
process A{E{t)) in (15. 5 p is given by 



(5.6) Q(x,t)= P{x,u)h{u,t)du 

Jo 

where h{u,t) is the density (15. 4p of the inverse subordinator E(t). The distribution 
function Q{x,t) solves the generalized Cauchy problem 

(5.7) i}D{dt)Q{x, t) = -^a{-iD,)Q{x, t) + H{x)(t)D{t, oo) 
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in the mild sense, where H{x) = I{x > 0) is the Heaviside function. Furthermore, 
P{x, u) solves the Cauchy problem 

(5.8) dtP{x, t) = -^^(-tD,)P(x, t)- P(x, 0) = H{x), 
and h{x, t) solves the inhomogeneous Cauchy problem 

(5.9) d:,h{x, t) = -ipD{dt)h{x, t) + S{x)(l)D{t, oo). 

Proof. The proof is similar to [311 Theorem 4.1]. Equation (15. 6p follows from a simple 
conditioning argument. Apply [3^1 Theorem 3.6] to see that Q{x,t) has FLT 

St f „-ifc.x^.„,_„ .N „,. 1 ^d(s) 



(5.10) Qik,s)= e-^' / e-''-^Qidx,t)dt ^ ,,, ^ ,, 

Jo V s ijAik) + iJois) 

and rearrange to get 

(5.11) iJD{s) Q{k, s) = -Mk) Q{k, s) + s-Vd(s). 
From [HH Eq. (3.12)] we get 

POO 

(5.12) / e-">z)(M,oo)rfM = s-VD(s). 



^0 

Now invert the FLT flSTTTD . using fl5TT2|) and J e-'''''H{dx) = 1, to arrive at flSTjl . 
It is well known that P{x,t) solves the Cauchy problem (15. 8p . see for example [20] . 
Equation (14.51) shows that the bivariate Laplace transform (LLT) 



oo /"OO 



-Xz-sttr- .^ ^. ^- 1 '^d{s) 



h(X,s)= / e-^'-''h(z,t)dtdz , , , ,. 

Jo Jo sA + ^b(s) 

This rearranges to 

Xh{^, s) = -ijD{s)h{\, s) + s"Vd(s). 

Inverting the LLT using (15.121) to see that h{x,t) solves (15. 9p . D 

For any random walk S{n) = X]r=i ^j' ^^^ compound Poisson process A{t) = 
S{Ni{t)) is a Levy process. Introduce IID waiting times (14.30 between these random 
walk jumps to get a CTRW. In this case, the CTRW is exactly of the form A{E{t)), 
without passing to the limit. Then the governing equations in Theorem 15.11 pertain 
to the CTRW itself. 

Theorem 5.2. Assume D{t) is a subordinator without drift such that conditions 
(14. 2 p and (15. 3p hold, and let E{t) be the inverse subordinator (II. 7p . Take J„ IID 
waiting times according to (14. 3p . and let Nij{t) denote the renewal process (14.71) . 
Take Y^ IID jumps on M.'^, independent from (Jn), with common distribution fi, and 
let S{n) = XliLi^i- Then the distribution function P{x,t) = ¥{X(t) < x) of the 
CTRW X{t) = S{N£,(t)) solves the generalized Cauchy problem 

(5.13) '^Didt)P{x, t) = -XP{x, t) + \ I P{x- y, t) fi{dy) + H{x)(j)Dit, oo) 

15 



in the mild sense. Furthermore, X[t) = A{E{t)), where A{t) = S{Ni{t)) is a com- 
pound Poisson process. 

Proof. Theorem O yields Noit) = Ni{E{t)), and then the CTRW is 

Xit) = SiNoit)) = S{N^{E{t))) = A{E{t)). 

A standard conditioning argument shows that the compound Poisson FT P{k,t) = 
g-t'/'A(fc)^ where the Fourier symbol ipAik) = \{l—jl{k)). The inverse FT oiijJAik)f{k) 
is 



(5.14) M-^D,)f{x) = -A/(x) + XJ f{x-y) fi{dy) 

using the FT convolution property. Now Theorem 15.11 implies that (15.131) holds. D 



Remark 5.3. In the situation of Theorem 15. 2[ where A{t) is compound Poisson, the 
distribution function P[x,t) = F{A{t) < x) solves the Cauchy problem (15.81) . which 
can be written in this case as 

/oo 
P{x-y,t)fi{dy); P{x,0) = H{x). 
'OO 

This is the Kolmogorov forward equation for the Markov process A{t). If /i has density 
w{x), apply dx on both sides of (15.151) to see that the probability density p{x,t) = 
dxP{x,t) of A{t) solves (II. 6p . If D is the stable subordinator with Laplace symbol 
^^(s) = s^, then fl57[3|) holds with ^^(^,00) = t-^/r(l - /3) and ^0(^4) = ©f, the 
Riemann-Liouville fractional derivative. The Riemann-Liouville fractional derivative 
is defined for < n — 1 < /3 < n by 

(5.16) Df^(t) = Y^^^^^^ I (t - rr-'-^g'^-\r) dr, 

which differs from the Caputo derivative (ll.5p in that the derivative is applied after 
the integration. The LT of I])^g{t) is s^g{s). Apply d^ to both sides of (I5.13P in this 
case to get 

f t~^ 

Dfp(x, t) = -Xp{x, t) + X p{x- y, t) fi{dy) + (5(x) ^ _ 

the fractional kinetic equation of Zaslavsky j49]. To recover (II. 4p . use d^g{t) = 
^t9it) - 9{0)t-^/T{l - f3) and p{x, 0) = 6{x). 

Remark 5.4. In the special case where /i = ei is a point mass, so that Yn = 1 almost 
surely, A(t) = Ni{t) is a Poisson process with rate A > 0. Then the distribution 
function P{x,t) of the renewal process Nuit) = A{E(t)) solves 

(5.17) Mdt)Pix, t) = -A[P(x, t) - P{x - 1, t)] + H{x)(t)D{t, 00). 
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If D is the stable subordinator with Laplace symbol ^d(s) = s'^, Equation f l5.17p 
reduces to 

d'^P{x, t) = -X[P{x, t) - P{x - 1, t)] 

as in Remark 15.31 The probability mass function p{n,t) = P{n, 1) — P{n — l,t) = 
AP{n,t) for n > 0. Apply the difference operator A on both sides to obtain 

dtP{n,t) = -A[p(n,t) - p{n - l,t)] 

as in Jumarie 



Remark 5.5. Scher and Lax |35] showed that a CTRW with waiting time distribution 
u and jump distribution z/ has FLT 



si— u{s)i){k) 



where i){k) = J e ^^'^v{dx). To reconcile with Theorem 15. 2[ recall from (14. 6 p that the 
waiting times (14. 3 p in Theorem 15.21 have LT 

u{s) = f e-''uj{dt) ^ 



X + ipois) 



and then it follows that ipnis) = A(l — u{s))/oJ{s). The jumps Yn in Theorem 15.21 
have Fourier symbol ipA{k) = A(l — fi{k)) and then (I5.10p implies 

Qik,s)-^ ^i.(.) _1 ^ _1 l-^{s) 



s^A{k) + i^D{s) s 1^ + (I - f,(k)) sl-u{s)fi{k) 

which provides a different proof that the CTRW equals A{E{t)) in this case. To 
simulate the sample paths of the non-Markovian process A{E{t)), it is sufficient to 
simulate the CTRW. In particular, the renewal process ND{t) gives the exact jump 
times of the inverse subordinator E{t). 

Remark 5.6. In the general case, where A{t) is not compound Poisson, Theorem 15.21 
provides a useful approximation. Given a Levy process A{t), take Yn = A{n) — 
A{n — 1), so that S{n) = A{n). Take N{t) a Poisson process with rate 1, so that 
5'(A~^A^(At)) is compound Poisson with Fourier symbol 

X{l-e-^^'^^^''^)^^A{k), asA^oo. 

Then S{X-^N{Xt)) => A{t) as A -> oo, and the CTRW with IID waiting times (gSD 
and these compound Poisson jumps converges to A{E{t)) as X ^ oc. As in Remark 
15.51 this fact can be used to simulate sample paths of the process A{E{t)). This fact 
has been exploited by Fulger, Scalas and Germano [TSj to develop fast simulation 
methods for space-time fractional diffusion equations. 
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Example 5.7. Tempered stable subordinators are theoretically interesting |11I13] and 
practically useful [IHIEZ]. Take D{t) tempered stable with Laplace symbol iPd{s) = 
(s + a)^ — a^ for a > and < /3 < 1, and let E(t) be its inverse (11.71) . Theorem 14.11 
shows that Ni{E{t)) is a renewal process. Let (r„) denote the arrival times of this 
renewal process, and use (14.41) to get 

^(^ "') = x + {s + a)^-aP- 

This tempered fractional Poisson process Ni{E{t)) has tempered Mittag-Leffler wait- 
ing times, but with a different rate parameter: Use (12. 4 p to see that the Mittag-Leffler 
PDF f{t) = dt[l — Ei3{—'r]t^)] has Laplace transform 77/(77 + s^), and so 

-st f- /j.\ „— at Tj. / 



e-''f(t)e-'''dt = - 
V + {s + a) 



Of course f{t)e "* is not a PDF, and in fact we have (set s = above) 

f{t)e-'''dt 








at/ 



Then the tempered Mittag-Leffler PDF fa{t) = f{t)e-''\ri + a^)/?/ has LT 



^0 r] + {s + aY \ + {s + aY-aP 

when Tj + a^ = X. Cartea and Del-Castillo [12] show that the tempered fractional 
derivative ipD{dt)g{t) = e~°* 9f [e"* g(t)]—a^g(t). It is also known (e.g., see [1]) that the 
corresponding Levy measure is exponentially tempered: ipoidt) = e~°''^ip{dt), where 
ip(t, 00) = t~^/r(l — /3) is the Levy measure of the standard /3-stable subordinator. 
Then Theorem 15.21 shows that the CTRW with tempered Mittag-Leffler waiting times 
and compound Poisson jumps solves a tempered fractional Cauchy problem 

g-at Q^^^at p^^^ ^^] _ a/3p(a;, f) = 4) A^-l D ^) P {x , t) + H{x)(j)D{t, Oo) 

with ipAi-iD^) given by (ISAij) and 0D(t,oo) = (3 J^°° e-^H-^-^dt/Til - P). More 
generally. Theorem 15 . 1 1 shows that the distribution function of the CTRW scaling limit 
A{E{t)) is governed by this equation, with the corresponding operator ipA^—iD^). 
Apply dx on both sides of (15.171) to see that the PDF of the renewal process with 
tempered Mittag-Leffler waiting times solves 

e-^^d^le^^'pix^t)] - a^p{x,t) = -X[p{x,t) - p{x - 1, t)] + 6{x)(f)D{t, oo). 

A wide variety of tempered stable models in M'^ are discussed in Rosihski [13] . Random 
walks in M'^ with tempered stable scaling limit are developed in [13]. For exponentially 
tempered stable waiting times in M^, a renewal process with tempered Mittag-Leffler 
waiting times gives the same process exactly, without taking limits. This can be 
useful for simulating sample paths. 
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Example 5.8. Chechkin et al. [13 [S] used distributed order fractional derivatives 
to model multi-scale anomalous subdiffusion, where a different power law pertains at 
short and long time scales, and ultraslow diffusion, for a plume of particles spreading 
at a logarithmic rate. Given a finite Borel measure v on (0, 1), the distributed order 
fractional derivative is defined by 

(5.18) W,g{t) = f d^gitMdP), 

Jo 

where 9f is the Caputo fractional derivative (11. 5p . If u is discrete, this is a linear 
combination of fractional derivatives. Let D{t) be the distributed order stable sub- 
ordinator with Laplace symbol iPd{s) = J s'^h'{d(3) and E(t) its inverse (II. 7p . Let 
i/(d/3) = piP)dp for some p G C^O, 1), then by (2.19) in Kochubei [21] 



(5.19) P{Jn >t) = E(e^^W) = - / r-ie"*"$(r, 1 

TT Jo 



)dr 



where 
Mr. 11 



/oV/^sin(/37r)r(l-/3)p(/3)rf/3 



[Jl rP cos(/37r)r(l - /3)p(/3)d/3 + A]^ + [/^ r^ sin(/37r)r(l - /3)p(/3)rf/3]^ 



Substitute (15.191) into (14. 9 p to obtain an explicit formula for the probability mass 
function of the distributed order Poisson process. 

If v{df5) = p{(5)dp, where p(/3) is regularly varying at /3 = with index a — 
1 for some a > 0, then ipnis) = R{logs) and R is regularly varying at infinity 
with index —a, see [35| Lemma 3.1]. Then E{t) is "ultraslow" in that E,{E{t)'') = 
S{\ogt), where S varies regularly with index 7a at infinity, by [351 Theorem 3.9]. 
Take an IID sequence of mixing variables (Bi) with distribution /i concentrated on 
(0,1), and assume '^{Ji > u\Bi = (3) = c~^u~^ for u > c~^'^, so that the waiting 
times are conditionally Pareto. Then [35l Theorem 3.4] implies that the distributed 
order stable subordinator is a random walk limit YlfJi ^i ^ -^(^)- This requires 
/(I — f3)~^fi{df3) < 00 so that h'{d(3) = r(l — f3)fi{df3) is a finite measure. An easy 
computation shows that the Levy measure 4>D(t, 00) = f^ t~^z/((i/3)/r(l — /3). Then 
Theorem 15.11 implies that a CTRW with these conditionally Pareto waiting times 
has a scaling limit A{E{t)) whose distribution Q{x,t) solves the distributed-order 
fractional diffusion equation 

3'',Q{x,t) = -^AHD^)Q{x,t). 

If A{t) is compound Poisson, Theorem l5. 21 shows that the distribution function P{x, t) 
of a CTRW with waiting times (14.31) solves 

B';P{x, t) = -\P{x, t) + X P{x- y, t) fi{dy), 
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without passing to the hmit. Then the PDF p{x, t) of the renewal process with waiting 
times ( I4.3P solves 

B'^p{x,t) = -X\p{x,t)-p{x-l,t)]. 
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